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Abstract 

We continue with the study of the Hankel determinant, 

D n (t,a,(3 ) := det x^ +k w{x\ t, a, (3)dx 

generated by a Pollaczek-Jacobi type weight, 

w(x\ t, a, 0) := x a (l — x)^e~ t ^ x , x G [0,1], a > 0, /3 > 0, t > 0. 

This reduces to the “pure” Jacobi weight at t = 0. We may take a € M, in the sit¬ 
uation while t is strictly greater than 0. It was shown in Chen and Dai (2010), that 
the logarithmic derivative of this Hankel determinant satisfies a Jimbo-Miwa-Okamoto 
a -form of Painleve V ( Py ). In fact the logarithmic of the Hankel determinant has an 
integral representation in terms of a particular Py. 

In this paper, we show that, under a double scaling, where n the dimension of the 
Hankel matrix tends to oo , and t tends to 0 + , such that s := 2 n 2 t is finite, the dou¬ 
ble scaled Hankel determinant (effectively an operator determinant) has an integral 
representation in terms of a particular Pup. Expansions of the scaled Hankel deter¬ 
minant for small and large s are found. A further double scaling with a = — 2n + A, 
where n oo and t, tends to 0 + , such that s := nt is finite. In this situation the 
scaled Hankel determinant has an integral representation in terms of a particular Py, 
and its small and large s asymptotic expansions are also found. 

The reproducing kernel in terms of monic polynomials orthogonal with respect to the 
Pollaczek-Jacobi type weight, under the origin (or hard edge) scaling may be expressed 
in terms of the solutions of a second order linear ordinary differential equation (ODE). 
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With special choices of the parameters, the limiting (double scaled) kernel and the 
second order ODE degenerate to Bessel kernel and the Bessel differential equation, 
respectively. 

We also applied this method to polynomials orthogonal with respect to the perturbed 
Laguerre weight; w(x;t,a) := x a e~ x e -4 / 1 , 0 < x < oo, a > 0, t > 0. The scaled 
kernel at origin of this perturbed Laguerre ensemble has the same behavior with the 
above limiting kernel, although difference scaled schemes are adopted on these two 
kernels. 
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1 Introduction 

The determinant of the n x n Hankel matrix, 

/ w(x)x j+k dx 

J L 

has an equivalent representation as the multiple integral [39] , 



where w(x) is a positive weight function supported on L (cK) and v(x) is known as the 
external potential. Hankel determinant a fundamental object unitary random matrix theory 
[39] with many applications, in mathematics, physics and other areas, for example, in wireless 
communications [I]. For instance, such determinant that arises from the singularly perturbed 
Laguerre weight, becomes the moment generating function of certain linear statistics ra- 
Such determinant also appears in the computation of the Wigner delay time distribution in 
chaotic cavities, studied from the point of view of large derivations[35]. We refer the Reader 
to mm for related material. 

Here are some well-known facts on orthogonal polynomials. 

The joint probability density function of the eigenvalues x\,x 2 , ■ ■ ■ ,x n of a nxn Hermitian 
matrix ensemble can found, for example, in [39] . 

11 n 

P(x i,x 2 ,...,x n ) = — Yl {xj-Xkf Y[w(x £ ), 

Un[w \ n - i< J<k < n i=i 


^ 0<j,k<n—l 


From which, he m -point correlation function, follows, see [23, [59] ; 

X 2 , • • • , 2-m) ' 

An equivalent expression reads, 

Rm{x\, ■ ■ ■, x m ) = det (. K n (x k , . 

Here kernel K n (x,y ) is defined by the monic polynomials P n (x) orthogonal with respect to 
the weight w(x) on L, 


ni 


n 


771 


p(x i, x 2 , ■ ■ ■, x n )dx m+ 1 • • • dx r 


' L n_m 


n —1 

K n {x,y) ■= V w ( x )V w (v)Yl 

3=0 


p j( x )Pji.y) 

h 3 


( 1 . 2 ) 


P n {x)P m {x)w{x)dx = h n 5 n 


(1.3) 
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and h n is the square of the L 2 norm. An immediate consequence of the orthogonality 
relations is the reproducing property; 


K n (x,y) = / K n (x,z)K n (z,y)dz. 


Further more, the monic orthogonal polynomials satisfy three terms recurrence relations, 


xPn(x) Pn+ 1(*^) ~t~ (X n P n (A) T (3 n P n —i (x), 


subjected to the initial data, Pq(x) = 1, and /3 0 P^i(x) = 0. 

With the aid of the Christoffel-Darboux formula [33], which is an immediate consequence of 
the three term recurrence relations, the kernel has a simple closed form 


tv- / \ / 7 T / 7 \Pn(x)P n — l(y) 

Kn{X,y) = V®Wv»(!/)-- 


- P n (y)P n _ i(g) 

(■ x-y) 


(1.4) 


It is of interest to investigate the feature of local eigenvalue correlation through a description 
of the correlation kernel, for n . For example, in the case of the Gaussian Unitary Ensemble 
(GUE), where w(x) = e~ x , x G R the limiting mean eigenvalue density reads, 


lim 

n—t oo y Tl 


R, i ^ \Z~2n: 


2 nx ) = 


2yT^ 


x A 


— 1<X<1. 


71 


This is the Wigner semi-circle law. Re-scaling with respect to a fix point Xq, leads to the 
sine kernel 


lim -JL, 

n ->°° V2n 


K n ( x o + 


nx 


--■> Xo + 


iry 


sin 7 t(x — y) 
n(x - y) 


yfin U yfhl, 

The limiting kernel becomes the Airy kernel with a suitable re-scaling at the edge \/2n, the 
edge of the eigenvalue spectrum, obtained by Tracy and Widom m, 


Kury(x, y) = lim -T—rK* ( v / 2n + -Ar. V2n + AT j - 
n ^°° 22/16 \ 22 n 6 22/16 / 


Ai(x)Ai'(y) — Ai\x)Ai(y) 


x-y 


Here Ai(z) is the Airy function. Tracy and Widom investigated the logarithmic derivatives of 
operator determinants, involving the Airy kernel, in the study of the level spacing distribution 
[46] , Chen and Ismail [l4j, obtained the Airy kernel by studying the limiting behavior of 
kernels generated in the situations, where v(x) is any even degree polynomials in x with 
positive coefficient in the highest order monomial. 

It is a useful technique to characterize the large n behavior of the scaled kernel in terms of 
differential equation. Tracy and widom mi adopted tools from integrable system to analyze 
the correlation kernel in the large n limit, and scale at hard edge of the Laguerre unitary 
ensemble (LUE). The limiting kernel can be expressed by the regular solutions of the Bessel 
differential equation. We restate the limiting kernel here, mi, 


A^Bessel (U 2/) 


<t>(x)y<j>’(y) - xcj)'(x)(j)(y) 
x-y 


(1.5) 
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where 4>{x) is the regular solutions of the Bessel differential equation ((2.14), [ 32 ]), 


x 2 (p"(x ) + x(f)'(x) + -(x — a 2 )(p(x) = 0 , ( 1 . 6 ) 

namely 0(x) = ^/JlJai^/x), g is a parameter, and J a (z ) is the Bessel function with order 
a, and 


A Bessel (-p £/) 


A Bessel (-G *^) 


A 

h 

4 


J*(Vx)VyJ' a (Vy) - VxJ'a(Vx)MVy) 

2(x - y) 


Ja+1 ( \/~x) Ja— 1 (V^)) ■ 


(1.7) 

( 1 . 8 ) 


This is known as the Bessel kernel. See [48] for further information. 

From the formulas of Laguerre and Hermite polynomial, Forrester |25j obtains the Bessel 
kernel and the Airy kernel after suitable re-scaling. Moreover, Nagao and Wadati [44] de¬ 
duced the Bessel kernel by scaling the Jacobi ensemble at the hard edges, ±1. Kuijlaars 
and Zhang [38] obtain a limiting kernel as a generalization of Bessel kernel by scaling the 
correlation kernel of complex Ginibre random matrices at the hard edge, see the references 
therein for more information. 

The theory of integrable kernels was put forward in [33]. A condition for a kernel K(yx,y) 
to be integrable, is that it can be expressed as the sum of functions fk(x) and gk(x ), that 
is, 

1 p p 

K(x,y) = - ^ fk(x)g k (y): where Y ^ f k (x)g k (x) = 0. 

x — y ' 

y k =0 k =0 

The Sine, Airy and Bessel kernel are all integrable. 

In our approach, we study the Hankel determinant directly, without expressing it in the form 
of det(/ n - K n ). 

For the problem at hand, the Hankel determinant reads, 


D n (t,a,(3 ) := det (^J x^ +k w{x\ t, a, (3)dx^J 


n —1 


j,k =0 


where 

w(x;t,a, (3) = x a ( 1—x) / 3 e _ ^ x , 16 [0,1], t > 0, /3 > 0, a > 0, (1.9) 

is the Pollaczek-Jocobi type weight. 

For t > 0, -y- 0, as x —* 0, with far greater speed than x a tends to 0, if a > 0. 

The same can be said, for a < 0, as long as t > 0. The Pollaczek-Jacobi type weight 
violates the Szego condition (see [43] and [TO]), which reads, 


r 1 |lnu;(a;)| 
/o a/x(1 - x) 


dx < oo. 
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Any monic polynomial orthogonal with respect to some weight can be represented by Heine’s 
multiple-integral, and in our case, 


P n (z]t,a,(3) 


1 

n\D n (t,a,/3) 



) n 

l<7<fe<ri 


n 

Xk) 2 n w(xi ; t, a, /3)dxg. 

t= l 


The constant term of our orthogonal polynomial has the closed form expression, 


(~l) n P n (0-t,a,(3) 


D n (t, a + 1, P) 
D n (t,a,/3) 


The remainder of this paper is organized as follow. In Section 2.1, our interest lies in 
a double scaling analysis, where t —> 0 + , and n —> oo, such that s = 2n 2 t is finite, 
we shall see later that that double-scaled and in some sense infinite dimensional Hankel 
determinant has an integral representation in terms of a particular P m '. Its logarithmic 
derivative satisfies a particular Jimbo-Miwa-Okamoto a -form Painleve. This double scaling 
analysis is based on results of Chen and Dai [TO] , obtained in the finite n situation. In 
Section 2.2, we obtained the asymptotic expansions of scaled Hankel determinant as formal 
series for small and large s . Furthermore, the constant in the large s expansion is found 
later, in Section 2.4. In Section 2.3, we introduce a new double scaling scheme, where 
t —» 0 + , n —> oo, p = j3n —* oo, with fixed P and s = 2(1 + P)n 2 t is finite. In this case 
the infinite dimensional determinant may be characterized by the same Painleve equations 
that appeared theorem 3, after a minor change of variable. In Section 2.4, an evaluation 
is made on the constant term of the monic orthogonal polynomials P„(0; t, a, P ), for large 
n , and s = 2 n 2 t. This comes from an application of the Szego limit theorem for Toeplitz 
determinants, but adapted to Hankel determinants. From these results, the constant c(a) 
appears in the asymptotic expansion of the double-scaled Hankel determinant, for large s 
is found. In Section 2.5, we combine Normand’s formulas |4D] with method of Forrester 
and Witte HZ], to find asymptotic expansions of the Hankel determinant for small s, in 
agreement with our expansion in Section 2.2. 

In Section 3, we propose another double scaling scenario, namely, a = —2n + A, 
t — > 0 + , n —>■ oo, where A < 0. In such a way that s = nt > 0 is finite. We note that 
our Hankel determinant under this double scaling scheme can be reproduced by another 
perturbed Laguerre weight, namely, 


w(x] t ) 


(x + t) ( P +/3 )x /3 e x *, A < 0, P > 0, t > 0, x G (0, oo), 


studied in [2D]. In this case, the Hankel determinant has an integral representation in terms 
of a particular Py , (equivalent to a Pm )• From which we determine its small and large s 
expansion. Moreover, the constant which appear in the asymptotic expansion of the double- 
scaled Hankel determinant, for large s, is found. 

In Section 4, combining the ladder operator relations in x obtained in [10] and further 
ladder operator relations in t obtained here, satisfied by the monic orthogonal polynomials 
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Pn (x]t,a,/3) (see theorem 12), we construct a Lax pair, involving derivative in x and 
derivative in t. The natural compatibility condition, reproduces certain results of [TQ] . In 
order to analyse the limiting behavior of the kernel arising form the Pollaczek-Jacobi type 
weight, a scaling scheme is introduced. Here, t —* 0 + , n —* oo, such that s = 2 n 2 t is 
finite. The “coordinates”, x and y have been re-scaled to x = y = ■ Ultimately, 

this shows that the limiting kernel may be characterized by solutions a second order ODE. 
If s = 0, this essentially reduces to the Jacobi weight, and the limit kernel and the second 
order ODE reduce to the Bessel kernel and Bessel differential equation, respectively. 

In Section 5, we adopt the method in Section 4 to study the kernel arising form the singularly 
perturbed Laguerre weight. We adopt another scaling scheme, where t —* 0 + , n —> oo, such 
that s = (2 n + 1 + a)t is finite. The “coordinates”, x and y have been re-scaled to x = j- 
and y = £■ 


2 Double scaling analysis I. 


Chen and Dai uni applied the ladder operator method to investigate the Hankel determinant 
obtained from the Pollacaek-Jacobi type weight. It was found that the logarithmic derivative 
of the Hankel determinant satisfies a particular Jimbo-Miwa-Okamoto a -form of Painleve 
equation. 

An immediate consequence of the relationships obtained in jT0], is that the Hankel determi¬ 
nant has an integral representation in terms of a Py transcendent in the variable y(t, a, (3). 
See Lemma 1. We shall be concerned with the behavior of the Hankel determinant, as n , 
the dimension of the Hankel matrix tends to infinity. For this purpose, a double scaling 
scheme is introduced, namely, sending n — > oo, t — » 0 + , and such that s := 2 n 2 t remain 
fixed. We recall theorem 5.4 in in- 

Theorem 1. The logarithmic derivative of the Hankel determinant with respect to t, 

H n (t, a, f3) := t^-ln ^ ,a ’^ = (2n + a + P)(r*(t) -r n (t)), (2.10) 

CLt L>n\U, o, p) 

satisfies the following ordinary differential equation: 

0 tny = [n(n + a + P)-H n + (a + t)H' n ] 2 + 4 H' n (tH' n - H n ){P - <), (2.11) 


with the initial data H n (0,a,/3) = 0. 

Here r*(t) and r n {t) variables defined in [TO] , 

By a change of variable, the above ODE turns out to be a particular Jimbo-Miwa-Okamoto 
a -from of Py . See in for more details. 

We now recall the Theorem 7.2 in [10] and replace S n (t) in that paper by y(t,a,(3). 


Theorem 2. Let 


y(t,a,/3 ) 


Rn(t) 

2jTi -\-1 -\- ex -\- /3 


( 2 . 12 ) 
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Then y(t,a,/3) satisfies the following differential equation: 


V 


3y-l 2 y' (2n + l + a + ff(y - T) 2 y 
2 y(y - 1 r y) t 2t 2 


(y - x ) 2 /3 2 ay _ y(y + 1) 

2t 2 y + t 2(y-ly 

(2.13) 


which is a Py ((2n + 1 + a + /3) 2 /2, —(3 2 / 2, a, —1/2). T7ie boundary condition is r/(0, a, /3) = 

1. 


From results obtained in [10], we show that the Hankel determinant has an integral 
representation in terms of Py. See the Lemma below. 

Lemma 1 . The logarithmic derivative of the Hankel determinant H n (t,a, ft), and r* , can 
be expressed in terms of y{t,a,/3) and y'{t,a,ft) as follow, 


H n (t, ex, ft) 


t d_ ln D n {t,a,ft) 
dt 11 D n (0, a, ft) 


1 

*y{y- l) 2 


ft 2 + 2 y 3 - 2(2n + a + p) 2 y 3 


2 (t + P)(2n + a + ft)y 3 - y 4 - t 2 y' 2 


+ (2 n + a + ft) 2 y 4 — 4 n(t + /3)y — 2/3 (t + a + 2 /3)y — 2 ty'y + 2 ty'y 2 + 4n 2 ?/ 2 

+l/ 2 ((t + a) 2 + 4t/3 + 6 aft + Qft 2 — l) + 4n(2t + a + 3 ft)y 2 \ + n{n + a — t), 

(2.14) 


with the initial data H n (0,a,ft) = 0. 

Proof. Recalling identities (5.8), (5.9) and (7.2) in [TO] as, 

R (f) = (2n + 1 + a + ft) [2r 2 + (t + 2/3 - 2-r;)r n + (2n + «)< - nt - tr' n ] 
2[(r* - r n ) 2 + (2 n + a - t)r* + (/3 + t)r n - nt] 


1 _ 2r 2 + (t + 2/3 - 2r*)r n + (2n + a)r* n - nt + tr^ 

-Rn(t) 2(2n + 1 + a + /3)(/3 + r n )r n 

and 

< = ^ - (2n + 1 + a + P){2r n - R n + /3)] + r n - R " ^ 

Zlh n Z 

With the aid of above equations and the definition of y(t, a , /3) in (12.121) . the variables r n (t) 
and r*(t) may be expressed in terms of y and y'. Straightforward computation produces 

mm- □ 

Hence the Hankel determinant has an integral representation in terms of y and y' , and 
that y(t,a,j 3) satisfies the Py given by (12.131) . 




















2.1 Scaling limit of the Hankel determinant in terms of Painleve 
equations. 

Carrying out the double scaling and combining with Theorem 1, Theorem 2 and Lemma 
1, we find that the (effectively) infinite dimensional Hankel determinant has an integral 
representation in terms of (I 2.16 1) . Moreover, the logarithmic derivatives of such a Hankel 
determinant satisfies another a -form of the corresponding Painleve equation. 

Theorem 3. Let 


y(t, a , /3) := 1 + ^ , and s := 2n 2 t. (2-15) 

n 2 

t —> 0 + and n —> oo , such that s G (0, oo). 

If 

g(s,a,P) := lim f (—,a,p) , 

7WOO \Zn z / 

then g(s,a,/3 ) satisfies 


1 g s s 2 2s 4g’ 


(2.16) 


with the initial data g(0, a, (3) = 0, g'(0, a, /3) = 2-, The equation Ii2.16\) is Pur (8, 2a, 0, —1) 

If 

U(s,a,P) := lim H n (—,a,/3) , 

n—>oo V Zn z J 

then 1-L(s,a,(3) satisfies, 


(s«") 2 + 4 (H') 2 UH' - H) - (an' + () =0, 


(2.17) 


with the initial conditions FL(t),a,P) = 0, TL'( 0,a,/3) = — ■ 

Furthermore, if 

A 7 a \ D n (s/2n 2 ,a,P ) 

A(s,a,^) := lim -, 

n^oo _D n (0,a,p) 

then 


n 


d 

ds 


In A(s, a, (3 ) 


(s S'' - #) 2 , 4as5f - s 


W 


+ 


16<? 2 



(2.18) 


Proof. Substituting the definition (12.15p into (1 2.13 1) . then we see that g(s,a,(3 ) satisfied 
Pm'(8, 2a, 0, — 1) . See (1 2.16 p . Plugging s = 2n 2 t into (1 2.11 p . we see that the limit of 
H n (t,a,f3) satisfies a a -from Painleve equation (1 2.171) . Moreover, substituting (1 2.15 1) 
into (1 2.14 1) . then the equation (1 2.18 1) is found. □ 
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Note that the differential equation is satisfied by g(s,a,/3) is a Pm'(8, 2a, 0, — 1), see 
EE302]. From P m /(8, 2a, 0, — 1) and the a -form of the Painleve equation, together with the 
boundary conditions, g(s,a,j3 ) and %(s,a,/3) are independent of j3 . We use these notations 
to distinguish this result from those previously obtained [9]. The cr -from Painleve equation 
of (I 2.17 1) has the same double scaling limit of the logarithmic derivative of the Hankel 
determinant generated by the singularly perturbed Laguerre weight, X a e~ x ~ t ^ x , where x > 
0, t > 0, and real a, however, the double scaling scheme is different from that in |9j. 

Remark 1: By the change variables 


F(x,a,P) = -g 

X 



1 


it is seen that the P m '(8, 2a, 0, —1) satsified by g becomes, 

+ ^ i 

a Pm(l, 2a, 0, —1). See [30]. The C potential (2.22) introduced in [9], after a change of 
variable, satisfies the above equation. 

There are three algebraic solutions of Pin'(8, 2a, 0, —1), 

1 2 

g(s,a,(3) = -s 3 , for a = 0. 


g(s,a,/3) 




for a = ±1. 


2.2 Asymptotic expansions of the scaled Hankel determinant. 


We assume that the solution of P ni '(8, 2a, 0, —1) for s —> 0 + has the power series ex¬ 
pansion ]Cylo a j s '\ with g(0) = 0, g'(0) = and substitute this into (I2.16L by some 
straightforward computations, one finds, 


g(s, a, (3) 


1 _ 1 2 3 

2a 2a 2 (a 2 — 1) 2a 3 (a 2 — 4)(a 2 — 1) 

5(—36 + lla 2 ) 5 

2a 5 (a 2 — l) 2 (a 2 — 4)(a 2 — 9)(a 2 — 16) 


3 9 — 6a 2 

' S + a 4 (a 2 - l) 2 (a 2 -4)(a 2 

+ (P(s 6 ), where a. j-7L. 


9) 


s 


4 


(2.19) 


For large and positive s , we assume that the solution of (12.161) has the following expansion 
3. The first term of the expansion is s 2 ^ 3 /2. A straight forward computation 

gives, 


1 2 

g(s,ot,P) =-s 3 


1) _i 
— S 3 + 


a 2 (a 2 


a 2 (a 


a i ala - 

— S 3 -j- - 

6 162 

2 — l)(2a 2 — 11) 4 

n - -s~3+0(s 


486 


1) 2 a(a 2 — 1) i 

— S"3 + V __ - 1 S _1 


486 


6561 


( 2 . 20 ) 
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Note that this solution becomes the algebraic solutions mentioned in the Remark 1, for 
a = 0 , and a = ±1 

In the next Theorem we obtain asymptotic expressions of the scaled Hankel determinant 
for small s and large s . 


Theorem 4. Under the double scaling scheme, the asymptotic expansions of the scaled Han¬ 
kel determinant generated by the Pollaczek-Jacobi type weight, has the following small and 
large s expansions: 

For small s, 


A(s, ct, f3) = exp 



10a 5 (a 2 — 


8a 2 (a 2 — 1) 6a 3 (o; 2 

(36 - lla 2 )s 5 
l) 2 (a 2 — 4)(a 2 — 9)(a 2 


— l)(a 2 — 4) 16a 4 (a 2 



3(2a 2 - 3)s 4 
— l) 2 (a 2 — 4)(a 2 — 9) 

( 2 . 21 ) 


where a ^ Z . 
For large s. 


A(s, a, (3) = exp 

a(l — a 2 ) 


9 2 3a i 1 — 6a 2 a(l — a 2 ) _x a 2 (l — a 2 ) _2 

c -S 3 -]-ss -|-InsH---—s 3 -]-i- -s s 

8 2 36 18 216 


486 


s + 0(s 3 ) 


( 2 . 22 ) 


where c = c(a) is an integration constant, independent of s . 
Proof. By (1 2.18 j) , we see that 


In A(s, a, j3) 



(&' - 9? 4 ajg - £ 2 

4 f,g 2 16 fg 2 


4<? + « 2 \ 

—) d( 


For small s , with g(s,a,/3 ) given by (1 2.19 11 the asymptotic expansion A (s,ck,/3), (1 2.21 1) 
follows immediately. Similarly, for large s, with g given by (1 2.20 11 , the equation (1 2.22 ]) 
is obtained, following straightforward computations. □ 


Note that the scaling limit of Hankel determinant via from the Pollaczek-Jacobi type 
weight is independent of P. 

In what follows, we give an account which will ultimately determine c(a). We start from, 

t 'i \ n t) frt. 4. _. m D n (t,a + l,P) 

( 1) Pn{0,ti a, p) — . 

Un\t'> P) 

Corollary 1. Sending n — > 00 , t — » 0 + and such that s = 2 n 2 t, is finite, then 


-l) n P n (0; 2^2 iOc,P) A(s, a + 1,/3) 


lim 

n^oo (—l) n P n (0; 0, a, P) A (s,a,P) 


= exp ci + -ss 


1 + 2a 
6 


In s 


a (a + 1) i 


6 


s 3 


0(0 T l)(2a T 1) _ 2 0(0 T 1) _ 


108 


* « - — s' 1 + O(s-s) ] , (2.23) 
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where c i = Ci(a) is a constant, independent of s, and 


Ci(a) = c(a + 1) — c(a ) 


(2.24) 


and c(a ) zs f/ze constant in \2.22\) . 


Proof. From the fact 


(- 1 ) 



(—l) n P n (0; 0, a, (3) 


,/d) D n (s/(2n 2 ),a,jd ) D n (0,a + l,/3) A (s,ce,/d) 


the equation (12.23p and (I2.24p follow from (12.22(1 . 


□ 


A computation that produces the constant ci(a) can be found in Section 2.4. 

2.3 The Hankel determinant for large ( 5 . 

In this subsection, we are interested in the behavior of the Hankel determinant for large fd 
and introduce a different scaling process. Let n —» oo, t —» 0 + , /3 := n/3, s := 2(1 + /3)n 2 t, 
such that s and f3 are fixed, we then obtain the same Painleve equations in the theorem 3; 
just replace f(t,a,/3 ) of (1 2.15 p in the theorem 3 with f(t,a,/3)/( 1 + j3). We state these 
results in the theorem below. 

Theorem 5. Let 


fd := nfd, s := 2(1 + (d)n 2 t , and y{t) := 1 + Q ’ ^ 



(2.25) 


t —» 0 + and n —> oo such that s and (d are finite, s £ (0, cxd) and (d £ (—1, oo). 

If 



then g(s,a,fd ) satisfies the following Pm'(8, 2a, 0, — 1), 



(2.26) 


with initial conditions g(0 , a, /?) = 0 , (?'(0, a, /3) = A-, 


If 



then T-L{s,a, fd) satisfies, 



(2.27) 
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with initial conditions "H(0,a,/3) = 0, Moreover, if 


D r 


A(s, a, j3) := lim 


2 ( 1 +/3)n 


a , n 


21^1 


D n ( 0, a, n, 


then 


Pi{s, a , /3) 


d_ 

ds 


In A(s, a, /3) 


(*/ ~ 9 ) 2 

V 


+ 


4ctsg — s 2 

16 p 



(2.28) 


Proof. Substituting 02.251) into the (1 2.13 p . g(s,a,/3 ) is found to satisfy (1 2.26 1) . Plugging 
[3 — ri/3 and s = 2(1 + /3)n 2 t into the a-from Painleve equation is satished by H n (t,a,j3) 
in the theorem 1, we see that 'H(s,a,/3) satisfies (1 2.27 p . Moreover, substituting (1 2.25 p 
into 0 2.14 p , the equation 02.28 p follows. □ 

Remark 2: Comparing theorem 3 with theorem 5, one finds that the Painleve equations 
are the same, although we emphasize that their scaling scheme are different from each other. 


2.4 Large n behavior. 

In order to find the constant ci(a) in (12.231) . we need to determine the large n behavior of 
the constant terms of the monic orthogonal polynomial, namely, P n (0; t, a, /3). The Szego 
limit theorem which computes the determinants of the finite section of Toeplitz matrix with 
nice symbols, can be adapted to orthogonal polynomials on the line. 

The large n computation for the orthogonal polynomials where the potential v satisfies the 
convexity condition pa can be found in EnanE] and also in [34]. 

In fact, as n —> oo , P n {z) is approximated by 


P n (z) ~ exp [-Si(z) - 5 2 (2f)], 


(2.29) 


valid for z [a, b\. Here ^(^ and S 2 (z) are given by ((4.6) and (4.7) in [T9]). These 
formulas are 


exp(—5i(*)) = 


z — b 


z — a 


+ 


z — a 
z — b 


, z £ [a,b ]. 


(2.30) 


and 


So(z) = —nln 


y/z — a + y/ Z — b 


f v(i) 

y/(z-a){z-b) , 1 

Ja y/(b — x)(x — a) 

x — z 


dx, z ^ [a, b], (2.31) 
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For the problem at hand, the Pollaczek-Jacobi type weight of (11.91) . then v(x) and v'(x) 
are given by, 

v(x) = — lnrcfa;) =- oAnx — /31n(l — x ), and v'(x) = -- — -——. 

x x 2 x x — 1 

Substituting v(x) and v'(x) into the following identities, 

f b V(x) 


and 


{b — x) (x — a) 
xv'(x) 


dx = 0, 


y/ (b — x){x — a) 


dx = 27m. 


(2.32) 


(2.33) 


These can be found, for instance, in mm- From (I2.32p . (I2.33f) and integral formulas in 
the Appendix A, it is found that a and b satisfy the following algebraic equations, 


a A b a 
-5 ~t H- 


13 


and 


2 ( 06)2 \fab -^/(l — a)(l — b) 

t (3 


= 0, 


yfah a/( 1 — a)(l — b ) 


T 2n A ol T (3 — 0. 


(2.34) 


(2.35) 


Here the parameters a and b determines the end points of the support of the equilibrium 
density. 

Let X := l/\Zab, and eliminating a + b from (12.341) and (12.351) . then X satisfies the 
quint ic, 

XH f3 2 XH v Xt , n n 

-— — 7 - 77777 A cxX — — — (2 n A 01 A /3) — 0. 

2 2{tX A 2n A a A f3) 2 2 y 

We now state a theorem which describes the large n asymptotic of P n (0; t, ct, f3) , without 
displaying the detail steps involved, since these are quite straightforward. 

Theorem 6. If v(x) = — In w{x) = y — aln x — /51n(l — x), x G [0,1], t > 0, a > 0, 
(3 > 0, the evaluation at z = 0 of Si(z ; t, a, (3), S 2 (z] t, a, (3), and P n (z ; t, a, (3) given by 


exp [—Si(0; t, a, (3)\ ~ 2 l (2n A a A (3 ) 2 {2 1 t(2n A a A (3) 2 } e, 


(2.36) 


exp [—52(0; t, a, (3)] ~(—l) n 4 n (2n A 01 A /3)“exp {2 l t(2n A a A (3) — (2a A /3)ln2 


I 1 rv 

+ -{2~ 1 t{2n A a A (3) 2 }z - — In (2 _1 t(2?r A a A (3 ) 2 ) 


(2.37) 
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and 


P n { 0; t, a, (3 ) ~ exp [-^(O; t, a, (3) - S 2 { 0; t, a, f3)} 

"3 


(—l) n 4 _n n" + 5 2“^ + " + ^ 


exp 


a 


6 


In {2 1 t(2n + a + (3) } 


{2 t(2n + ol + /3) } 3 — — In {2 t(2n + cr + /3)" } 

Z O 

(2.38) 


T/ie a&ove asymptotic estimations are uniform with respect to t G (0, t 0 ], 0 < t 0 < 
a > 0, /3 > 0, n —>■ oo srtc/i that n 2 t is fixed. 

In order to derive the constant ci(a) in fl 2.23 1) . we still need to obtain P n (0; 0, a,/3), 
which is the constant terms of monic polynomials orthogonal with respect to the ’’shifted” 
Jacobi weight w(x]0,a, /3) = x Q (l — x)^, x G [0,1]. This can be found from the monic 
polynomial orthogonal with respect to Jacobi weight w(x) = (1 — x) a (l + x)@, x G [-1,1]. 
Taking a result from mi. we find that, 

T(n + 1 + a)r(n + 1 + a + (3) (— l) n A~ n n a+ ^2~^ a+ ^ + ^ y/2^ 


-Fn(0; 0, a, f3) = (—1)’ 
P n (l]0,a,/3) = 


T(a + l)T(2n + 1 + a + (3) T(a + 1) 

T(n + 1 + /3)T(n + 1 + a + f3) £~ n nP + ?2~( a+ P + & y/zir 


T(/3 + l)T(2n + 1 + ol + /?) 
where use has been made of the asymptotic formula 


r(/3 + i) 


T(n + 1 + a) ~ v / 27m n+Q+2 e n as, n —> oo. 

Remark 3: To derive the constant C\(a) in 1 12.231) . we recall the asymptotic estimation 
of P n (0; t, a , (3) , f)2.38j) . as 

P n (0; t, a, f3) ~ exp [-Pi(0; t, a, (3 ) - S 2 (0; t, a, j3)\ 

T(a + 1) 


T(a + 1) 


y/2n 


-exp 


{2 L t(2n + a + /l) 2 } 3 


ry 1 

— — In {2 _1 f(2-n + a + /l) 2 } - - In {2 -1 t(2n + a + /l) 2 } 
3 o 


-exp 


(■-1)' n 4~ n n a+ h~' ^ a+/3+ i> V2vr 

T(a + 1) 

In {2 _1 f(2-n + a + /3) 2 } - ^ In {2 _1 t(2n + a + /l) 2 } 
3 o 


In r(Q / l 1) + \{2~H{2n + a + P) 2 }^ 

V27T 2 


(2.39) 


Hence, 


P n (0;t,a,/3) 

P n (0;0,a,/3) 


a 


exp -{2 1 t(2n + a + P) 2 } 3 — — In {2 1 t(2n + a + /?) } 


-In{2 1 t(2n + a +/3) 2 } + ci(a) ) , 
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from which ci(a) is found to be 


i,i r(a + 1) . 


y/2n 


Replacing 2 H{2n + a + (3) 2 by s , then (12. 2 Ml) follows. 


By (12.241) one finds 



giving 



which we recognize to be the Tracy-Widom constant that appeared in the Bessel kernel 
problem. 

Here G(z) is the Barnes-G function. 

2.5 An algorithm of the Hankel determinant for small s. 

In this subsection, we use the formulas derived by Normand in [40j and similar process in 
[57] to compute the Hankel determinant. We look for expansion of the Hankel determinant 
around t — 0. We scale the variable t as s = 2 n 2 t in the expansion and let n —> oo, in 
this way we verify the asymptotic expansion for small s derived in section 2. 


Recall that 


D n (t, a, P) = det [/^ +j (t)]^ 0 , 


where the moments Hi+j(t ) are given by 



Although the moments Hi+j(t ) depending on a and f3 , we do not display this to lighten 
notations. 

Lemma 2. The moments are, 


where (p m (t) and are analytic at t = 0, 

.. e~ t T(fd + l)T(m + a + 1) 



m e N, (2.40) 



iF^l + p, -m - a,t ), 


and, 


ipmid) = e < T (—m — a — 1) iFi(m + 2 + a + p, m + 2 + a, t), 
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where iF^a, b; z) is the confluent hypergeometric function of the first kind. 
Furthermore, n m (t) has series expansion around t = 0, 

Fm(t) = <Pm{0) + 0) + ^jVm(O) + Oft 3 ) + t" t+ " +A+1 '(/v(0) (1 + Oft 2 )) 


where 


Tm (0) 


£) ^2(m+a+A+l)^ 

r(/3 + l)r(m + a + 1 ) 


r(m + 2 + a + ft) 


tfn(0) = 


r(/3 + l)r(m + a) 
r(m + 1 + a + ft) ’ 


(2.41) 


(2.42) 


M o) = 


T(/3 + l)T(m + a — 1) 
T(m + a + /3) 


and 


(2.43) 


-0m(O) = T(-m - a - 1). 

The expansions are valid for a > 0, /3 > 0, m G N and | argi| < 7 r. 

Proof. From ([29], P367), we find 

. . f m n / , r —L 1 m+H . ( m 2 cx 2B m 1 -|- o \ 

Fm(t) = I x x (l-xfe *dx = t 2 e ^T(/3 + 1)VF (---,--- ,tj, 

where W denotes the Whittaker function 

r(-2ju) 


lF(A,/i, z) 


r(i-//-A) 


^+ 2 e -2 lFl (/r- A + -,2/r + l,z) 


T(2u) _ 1 1 _z 1 

—y—--Z 2 e 2 1 F 1 (— p, — A + — 2p + 1, z), 

1 U + T ~ A > 1 


subject to |argz| < 7r ([29], P1023). Hence fl 2.40 p and fl 2.41 p follow. 

Theorem 7. The series expansion of the Hankel determinant is given by, 

nfn + a + fl)t 


□ 


D n (t,a,P) = det [pk+j(t)}l/ =0 = D n (0,a,fl) 

t n(n + a +/3)(n(n +a +f3)a +(3)t 2 ^ 

+ 2a(a 2 - 1) + [ } 

T (n + a + l)r(n + 1 + a + fl)nt a+1 


a 


T(fl + l)T(a + l)r 2 (a + 2)sin(a7r)(n - 1)! JTLiC# + j) 


- (1 + oft)) + e> (t 2(a+1) ) 


(2.44) 
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subject to a ^ Z, a > 0 and | argt| < 7 r. .D ri (0, a,/3) has a closed form expression 


' Q+/3+1 \ /~y2 ( a+/3+l ^ ^ Qi-\-(3 ^ 


-p / a+£±± \ f £*+£ 

D fO a /H = 4 _n ( n+Q+ ^f27r')”—^_-_-_-_ - _ v 2 _ 

n[U,a ’ P) 1 j G(a + P + l)G(a + l)G(P + l) 

G(n + 1 )G{n + a + 1 )G{n + p + l)G(n + a + /3 + 1) 


x 


G 2 (n + !^±|±i)G 2 (n + 1 + ^-)T(n + s±|±l) 


, (2.45) 


where G(z) is the Barnes G -function. 
Proof. By ([230 , 


D n {t,a,f3 ) =det [/x fe+j (i)]^ 0 ~ det 


^fc+j(0) + ^fe+j(0) + t- 


■ Pk+j(°) 


2 ! 


+ ^ +Q+ V fc+j (o) 


n —1 


- fcj=o 


det bi+j(0)]tj, 0 + ( < ( > det [v3fc+j(0) + «vi +j (0)] 

, 2 rf+i(0) 


n—1 
k,j =0 


+ t 2 < t 2 > det 


^fc+j(0) + tp'h+ji®) + ^ 


2 ! 


n—1 


k,j =0 

n—1 


+ t“ +1 < t“ +1 > det [^ fc+j (0) + t k+j+a+1 ^ k+j (f))]l . 


(2.46) 


Here < t m > fit) denotes the coefficient of t m in the series expansion of f[t) in t. We 
note here Normand’s formula 


det 


r (zj + i) 


n —1 n —1 

=n 

*J=o j =o 


(6+(n-l- j)(l - a)) j r(^) 
T(aZj + b + n — 1) 


II ( z J~ z i)> 

0<i<j<N—l 


r (azj + b + i) 
which can be found in id. and 

det [ai, • ■ • , a j + bj, ■ ■ • a n ] = det [a 1 , • • • , a,, • • • , a n ] + det [ai, • • • , b jr - ■ , a„ ], (2.47) 

where a^ , b 7 are column vectors. We obtain, after an extensive computations, while bearing 
in mind 

A»(0, a,P) = det [Pk+j(0)} n k ~l 0 , 

the following 

< t > det [(p k+ j(0 ) + Vfc+j(°)]fcjio = < t > det 


r(^ + l)r(a + l+j + fc) 
T(2 + a + P + j + k) 


—t 


T(P + l)T(a + j + fc) 

r(i + a + p + j + /c) 


-i n—1 


k,j=0 


With mm , we find, 


<t> det 


T(P + l)r(a + 1 + j + k) r(/3 + l)r(a + j + k) 


T(2 + a + P + j + k) 


— t 


r(i + a + p + j + k) 


71—1 


k,j =0 


n(n + a + /3)t 

CK 


.D n (0,a,/3), 


The terms < t 2 > fit ) and < t a+l can be similarly derived, although with greater effort. □ 
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Corollary 2. From the series expansion of D n (t, a, j3) around t = 0 , see fl 2.44 1) , the series 
expansion of H n (t,a,/3) around t = 0 follows, 


TT ,, m ,d D n (t,a,(3) n(n + ot + /3) n(n + a + (3)(n(n + a + f3) + a/3) 2 
H n (t, a, fi) = t- In Dn(0jaii8) =--- 1 +- a>2 _ 1} -* 

+ o(f) + C+ 1 7 i’ ( " r 3 ? + 4 !n (n , + Mr" m d + °«) + ° (* !< “ +1) ) • 

(a + 1) 2 1 3 (a + l)sin(a7r)l (n)l [n + p) 

(2.48) 


subject to a ^ Z, a > 0, 

Proof. With H n {t,a,f3 ) 

(1 2.44 1) , and the identity 


t > 0 and | argi| < 7r. 

given by fl 2.10 1) . the series expansion of D n (t,a, j3) , equation 


T(n + f3) 

r(/3 + i) 


n —1 

+ j), 

3 = 1 


the equation (1 2.48 p is obtained. □ 

Corollary 3. Sending t —» 0 + , n —>■ oo, and s = 2n 2 t such that s E (0, oo) is finite, we 
have, 


n „ a\ v d , d S S 2 

H(S ’ /?) = h’i ln O„(0, a „g) = S S A(S ' ' 51 = "25 + - 1) 


7TS 


,a+1 


-(1 + C(s)) + C(s 2 (“ +1 )), 


+ C(s 3 ) 

(2.49) 


2" +1 (a + l) 2 T 3 (a + l)sin(a7r) 
subject to a ^ Z, a > 0 and | args| < 7r. 

Proof. Substituting s = 2n 2 f into (1 2.48 jf , and taking limit n —» oo, then (1 2.49 j) is 


obtained with 

q a+l 

lim 


r(n + a + l)r(n + 1 + a + (3)tt 


US' 


a+ 1 


n->oo 2“ +1 n 2 ( a + 1 ) (a + l) 2 T 3 (a + l)sin(a7r)r(n)r(n + /3) 2“ +1 (a + l) 2 T 3 (a + l)sin(a7r)’ 

where a ^ Z, a > 0 and T(n + a + 1) ~ T(n + l)n", as n —>■ oo, is used in the proof. □ 


Remark. Note that the small s expansion of "H(s,a,/3) is independent of j3 and 
coincides with the characterization of Painleve equation in the Theorem 3. In particular, 
the “first portion” of fl 2.49 1) is the same as the series expansion, 


njf m d , A 7 Q\ N' “ df , 4 (XSg - S 2 

PL{s, a, (3) = s— ln A(s, a , (3) = — 


+ 


1 

2 a 


16 g 2 
1 


9 ~ ~r 


cc 

4 


4a 2 (a 2 — 1) 2a 3 (a 2 — 4)(a 2 — 1) 


+ 0(s 4 


where the small expansion s of g(s,a,(3 ) is given by (1 2.19 p . We note that the “special 
function portion 1 ' of fl 2.49 p , cannot be obtained through power series. This method, pro¬ 
posed by Normand, capture both “special functions” and “power series”. T complexity in 
the computations increases very rapidly, as one attempt to include term with higher powers 
of s and should be a worthy future project to explore. 
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3 Double scaling analysis II. 

In this section, we propose another double scaling scheme by setting a = —2n + A, s — nt, 
as n —> oo, t —> 0 such that A and s are finite, the scaled determinant has an integral 
representation in terms of a particular Py, which can be reduced to a Pm, from which its 
expansions for small and large s can be found. 


3.1 The Painleve equations. 

For convenience, we introduce a number of items; 


d A n (t,a,P) 
dt 11 A n (0 ,a,P) 

where 


t d_ ln D n (t,a,P) 
dt n D n (0,a,/3) 


+ n 


n\, and H n (t, a, /3) := H n (t, a, (3) + n 2 — nX, 

(3.50) 


H n (t, a, f3) 


^d_ D n (t,a,f3 ) 
dt D n (0,a,/3) 


We proceed as follows. 

Substituting a = — 2 n + A into (13.5011 . we find, 


j.n 2 —n\ 

A n (t,A,/3) =- —D n (t,\,P) 

n! 


j.n 2 —nX 


n\ 



n 

JJ (Xj - x k ) 2 \Y x i 2n+X (1 - Xi) p e~~idx h 

1 <j<k<n i =1 


With the change of variable z = ——, mapping the interval [0,1] to [0, oo] , one finds, 

-l) n f 


D n (t, \,P) = 


n\ 


n ( z kZjf \\(zt + t) (A+/3) zfe t dz e , (3.51) 

(°,o°) n 1 < i<fc < n 


t=\ 


and 


(—l) n A n (0, A) 


G(n + l)G(n + 1 « A) 

GRj 


where A < 0. 

This leads us to study the Hankel determinant originates generated by the weight 

w(x; t) = (x + f) _< ' A+/3 ^ /3 e _x ” t , A < 0, j3 > 0, a: G (0, oo). 

We would like to consider a more general weight of the form, 

w(x; t) — (x + f)^ < - A+/ 3 ^x 7 e _:r_ *, A< 0 ,/ 3 > 0 , 7 > 0,7 — A — (3 > —1, x G (0, oo), (3.52) 

which the reduces to the one above (13.521) . if 7 = f3 and degenerates to the Laguerre weight 
at t — 0 . 
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We introduce two quantities in term of the monic polynomials P n (x) orthogonal with 
respect to 03.521) . over [0, oo). These are 


R n (t) := 


i_ r°°pM 

hn JO V 


w(y;t)dy , 


r n (t ) := 


7 r Pn-i(y)Pn(y) 


hn- 


n- 1 Jo 


y 


w(y,t)dy. 


(3.53) 

(3.54) 


Moreover, we find r n and R n satisfy the following Riccati equations, 


tj;r n = ( 2 n - A - p + 7 ) r n +n(n-X-/3)-{ 
dt R 


r n (r n - 7 ) r n (r n + 2n - A - P) + n(n - A - f3) 


1 - Rn 


— 2r n — 7 + (2n + 7 - /3-A + t(R n - 1)) R n . 


Let 


fr , n \ . d , D n (t, 7 , A,/9) 

#n(R 7 , A,/3) := t—In 


n—1 


*“'d„(0,7,A,/3) 


(3.55) 

(3.56) 

(3.57) 


where the derivation of — t YJj=o obtained through ladder operators, is not re-produce 

here. 

Note that, 


(—l) n A„(0, 7 , A, /3) = 
and a simply identity 


G(n + l)G(n + 1 + 7 — P — A) 


G(j — P — X) 

?n(t) = H' n (t). 


(3.58) 


Theorem 8. For finite n, let 


R n {t) 


1 

' 1 + y{t,i,\P) - 1’ 


t/ien 2 /( 7 , 7 , A, /3) satisfies 

^ (3y-l )(^) 2 F , (y-l ) 2 , ^ 2 - 7 2n \ (2n + l + 7 -/?-A) y y(y + i) 

9 = 2y(y— 1) -T + ^-r + A)!, --rj-(-i(Fi)’ 

with the initial conditions 


2 /( 0 , 7 , A, £) = y 1 (0, 7 , A, 0) 


7(2 n + 1+7 — P — A) 

(/3 + A)((/3 + A-7) 2 -l)' 


The quantity H n (t, 7 , A, /3) defined in (1 3.57 j) . and satisfies the following second order ODE, 
(tfitf = 4 (» + H' n ) (S„ ~ tH'f (H' n ~ 7 ) + ((« - 7 + /3 + - H„ - njf , (3.60) 
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with the initial conditions 


H n (0, 7 , A, P) = 0, H' n { 0 , 7 , A,/3) = 


77-7 


/3 + A — 7 ' 


Moreover, 


H n {t,l,\P) ■= 




m 


4y (y - l) s 


((/? + A)y - ^ + + A - 2 n)?/ 7 * 


£ 2 y 


4y 


2 (y — 1 ) 2(?/— 1 ) 4(y — l) 2 ' 

(3.61) 


Proof. Eliminating r n from (1 3.55 ]) and (1 3.56 p . substituting R n (t) = 1 + l/(y(t) — 1) into 
the resulting equation, we see that y(t) satisfies (1 3.59 1) . 

Eliminating r n from (1 3.58 1) and (1 3.55 Ih we obtain the following 

R n = F(H n ,H' n ,HZ), 

where F(-, ■, •) is a function of three variables, which we do not display. 

Combining the expression for R. n , (1 3.57 11 . (1 3.58 p , and 

_ R n ( n(n - A) + (7 - (5)n + tr n ) - r n (t + A + j3 - 7 - 2n) - 777 + r n (r n - 7 ) 


3 =0 


t( 1 - R n 


tRn( 1 - Rn) 

(3.62) 

we arrive at the non-linear second order ODE which is satisfied by H n (t, 7 , A,/?). 

Finally, combining (1 3.62 p and (1 3.57 p , we see that, 


H n (t, 7 , A, /3) = 


R n (77,(77, - A) + (7 - ff)n + tr n ) - r n (t + A + /3 - 7 - 2 n) - 777 f n (f n - 7 ) 


The equation (1 3.61 1) is found, with (1 3.56 p and R n (t) = 1 + l/(y(t) — 1). □ 

Remark 4: If H n {t , 7 , A, (3) —>• H n {t , 7 , A, (3) — nt — ?r( / d + A), then the non-linear second 
ODE (I3.60P becomes, 

fi/cV = (u n - tH' + 2 ( H') 2 - <uy 2 - 4 (tfl - 77) - p -\) (k - n - 7 V 


where = 277+ 7 + /? +A. The equation above is the Jimbo-Miwa-Okamoto [ 361 HI] a -form 

of Pv, with parameters 

V 0 = 0, Vi = — 77 , V '2 = -f3 - A, v 3 — -77 — 7 . 

Similar to the investigation in previous sections, let, 

R (s, 7 ) := lim H n (-, 7 , A, /?) , A(s, 7 ) = lim /n’ 7 ’ \ » ?( s > 7) : = lim V (-> 7 , A, /3 

n->oo \77 / n^oo A n (0,7,A,pJ n->oo \77 

(3.63) 


we have, 
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Theorem 9. Let s := nt , t —> 0 + , n oo, such that s is finite. Then g(s, 7) satisfies 
6q - l 5 Q (P -+■ A)“<71 Q ” i ' 2 9 ^ " ' 2 

.9 = 


2 ?(?-l) 

with the initial data 

£( 0 . 7 ) = 


(^) 2 _? + ^ +A ) 2 ^^- 1 ) 2 T /2 (^- 1 ) 2 


2 s 2 


2s 2 ? 


7 


?( 0 , 7 ) = 


27 


P + X y v ’ ,y (? + A)((? + A — 7) 2 — 1) 

The quantity TL(s, 7) satisfies a non-linear second order ODE, 

7/A 2 = 4 £/ fn-sil') + ((_ 7 + / 3 + A)H'- 7 >2 


sH! 


with the initial conditions 


£(0,7) = 0, #(0, 7 ) = 


7 


P + A - 7' 


Moreover, 


n = 


(sy*) 2 


{{P + *)y- if s v 


(3.64) 


(3.65) 


(3.66) 


4?(?-l) 4? y-1' 

The equation fjd.6Jf\) is satisfied by Py((/3 + A) 2 /2, — 7 2 / 2 , —2,0)). 

Proof. Inserting s = nt into (1 3.59 1) . (1 3.60 1) . (1 3.61 li and followed by the double scaling 
limit, the equations (1 3.64 1) . (1 3.65 11 and (1 3.66 II are found. □ 

Remark 5: If TL(s, 7 ) —> TL(s , 7 ) — s, then TL'{s , 7 ) —* TL'^s, 7 ) — 1 and the non-linear 
second order ODE of (I3.65j) becomes 

sH"\ 1 = ah’ (n - sH’) in' - 1 ) + f (— 7 + p + A) u' - p - a) 2 , 


This turns out to be a a -form of Pm . 

Remark 6: Substituting H n (t,\,(3 ) = H n (t, A, (3) + n 2 — nX and s = nt into the a -from 
Painleve equation (12.11 H . we see that TL(s, A, /?) satisfies another a -form of Pm, 

(sH") 2 = AH'(1 + n')(H - sH 1 ) + (A n' - fi) 2 . (3.67) 


3.2 The asymptotic Expansions of Pv((0 + A) 2 / 2 ,- 7 2 / 2 ,- 2 , 0 ) 

For small and positive s , we find, 


7 , _^7_ 

? + A (/3 + A)((/3 + A — 7 ) 2 — 1) 

{{P + -^) 2 ~ 1)(3(? + A) - 5q) + 7 2 (/3 + A + 7 ) 2 3 

(? + A) (/3 + A — 7 )((/3 + A — 7 ) 2 — l)((/3 + A — 7 ) 2 — 4) 1 1 
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where 0 + A — 7 0 Z. 

The solution of P v ((0 + A) 2 /2, — 7 2 / 2 , —2, 0) for large s . 
For large and positive s, we find, 


9(s) 


2 1 4 7 2 -1 

0 + A + + 16(0 + A) 
(4 7 2 - 1)(4 7 2 - 8(0 + A) 
256 


1 

2 + 
2 


1 - 4 7 2 , (4 7 2 

16 + 

^-s~ 2 + 0{s~^. 


l)(48(/3 + A ) 2 - 4 7 2 + 25) _3 
1024(/3 + A) S 2 


(3.69) 


Theorem 10. The scaled Hankel determinant generated by the Pollaczek-Jacobi type weigh 
has the following small and large s expansions: 

For small s, 


A (s, 7 ) = exp 


7 


-s + 


7(0 + A) 


0 + A — 7 ~ ' 2(0 + A - 7 ) 2 ((0 + A - 7 ) 2 - 1) 
_2 7 (0 + A)(0 + A + 7 )_ 

3(0 + A - 7 ) 3 ((0 + A - 7 ) 2 - 4)((0 + A - 7 ) 2 - 1) 


s 3 + £>(s 4 ) 


(3.70) 


where (3 + A — 7 ^ Z. 
For large s, 


A(s, 7 ) = exp ^c 2 - s - 2(0 + A)s 2 + 
(4 7 2 — 1)(0 + A ) 2 


(0 + A)(0 + A-2 7 ) ln: | (4 7 2 


64 


+ 0(s 


1)(0 + -M -i 

16 


(3.71) 


where c 2 = c 2 ( 7 ) is an integration constant independent of s 
Proof. From (1 3.57 p . (13.63 p and (13.66 p . we obtain. 

In A(s, 7 ) = 


r ( mt )) 2 

((0 + A)?(O-7) 2 

m) \ 

lo Umim-v 2 


1 

1 — 1 


df 

z ' 


Substituting (13.68P and (13.691) into the above formula, the equations (j3.70[) and (13.71 jl are 
found. □ 


At the end of this section, we evaluate the large n behavior of P n (0;f, 7 ). In fact, 

( 1 \n r> fn. 4 D n (t, 7 + 1 ) 

(-1) P„(0,t, 7 ) = ■ 

D n {t, 7 ) 

For the weight, 

w(x; 0) = x 7 ~( A+ ^e~ x , 7 — A — 0 > —1, x G (0, 00 ), 

the constant term has a closed form expression in terms of Gamma functions, 

(—l) n P n (0; 0, 7 ) = F(n |^ + 7 7 A ~ VTk n n+7 ~ a ~^ +l ^ 2 e _n , n 00 . (3.72) 

1 (1 + 7 - A - 0) 
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Corollary 4. With the double scaling, one finds, 


lim 

71—> OO 


(—l) n P n (0; t, 7 ) 
(—l) n P n (0; 0, 7 ) 


A(s ,7 + 1) ( t \ 0 + A (1 + 2^y)(/3 + A) _x 

-k(^) =exp l C3(7)+ ^ lns+ — t — - s 2 


(1 + 27 )(/3 + A ) 2 

16 S 


+ 0(s 


(3.73) 


where c 3 ( 7 ) a constant independent of s. Moreover, 


c 3 (7) = c 2 (7 + 1) - c 2 (7), (3.74) 

where 02 ( 7 ) is f/ie constant in (1 3.71 1) . 

Proof. From the fact that 

lim (~l) n Pn(0; 7 , 7 ) _ lim Ax(g/n,7 + l) D n { 0 , 7 ) _ A(a,7 + 1) 

n^-oo (—l) n P n (0; 0, 7 ) n ^oo D n (s/n, 7 ) Ai(0 ,7 + 1) A(s, 7 ) 

(13.73(1 and (13.74(1 follow. □ 


3.3 Large n behavior of P n (0;t,j). 


In this subsection, we obtain the constant, 03 ( 7 ), following similar steps in section 2.4. We 
substitute v(x,t) and v'(x,t) into (1 2.32 1) and (1 2.33 p . With the aid of the integral formulas 
in the Appendix A, we find a and b satisfy two algebraic equations, 


X + fi 

y/ {t + a) (t + b ) 



0 , 


(3.75) 


A + fi - 7 + 


a + b 
2 


t(A + /3) 

\/(t + o) (t + 6) 


2 n. 


(3.76) 


The parameters a and b , are, as before, the end points of the support of the equilibrium 
density. 

Let A := yfab , eliminating a + b from the above equations, we see that X satisfies a 
quintic equation, 


A 2 + 2t(2n + 7 - A -fi+^-t 2 ) (7 - A ) 2 - A 2 (A + fifi = 0. 


Note the large n of A" asymptotic of A", given by, 

7(0 +A), . 1 , 7(0 + A ) 2 , 7(0 + A )( 7 2 — 2({3 + A) 2 ) 


A ~ 7 


-(n£) 2 + 


■ (nt) + 


16 


nty 


(3.77) 


(3.78) 


with t G ( 0 ,p], 0 < ti < 00 , n —^ 00 such that nt is fixed. 

We now give the large n behavior of P n (0;f, 7 ), and place this in the Theorem below, 
without displaying the detail steps. 
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Theorem 11. If v(x) = — lnw(x; t) — (A + ( 3 ) ln(x + 1 ) — 7I11X + x + t, and 
A < 0,/3 > 0 ,7 > 0 , t > 0 ,7 — A — j 3 > — l,x € ( 0 , 00), the evaluation of the orthogonal 
polynomials at the x = 0, reads, 


■l)"-fn(0; t, 7) ~ n 


n+ 7 —a—/9+lg—ra 


exp 

(1 + 27)09 +A) 


lA , A + /?, , . 

7+ 2/ n7 + 7 + ~2~ 


(nt)-i - (1 + 27) ( /j + A) 2 (nt ) - 1 
16 


( 3 . 79 ) 


t/ie asymptotic estimation is uniform with respect to t e (0,11], 0 < + < 00, 
A<0,/3>0,7>0,7 — A — / 5 > — 1 , n —>■ 00 such that nt is fixed. 

Remark 7: To obtain the constant 03(7) in ( 13 . 731 ) . we rewrite the asymptotic estimation 
of (—l) n P n ( 0 ; t, 7), the equation ( 13 . 79 ft . as 

(—l) n P n ( 0 ; t, 7) ~ n n+7 - Q -^+ 5 e- n exp 

(I + 27XI +A) , _ (l+2 7 )(/? + A) 2 A 

4 v ' i 6 y 


1 A . A + j 3 . . 

7+ 2/ ln7 + 7 + ~^ 2 ~ ln (nt) 


Taking into account the extreme right of ( 13 . 721 ) . we see that, 


-Fn( 0 ;i, 7) ( A + [3 . . (l + 27)(/? + A) 1 (1 + 27) (/5 + A) 2 x 


exp c 3 + 


ln(?7i) + 


(nt)’ 


16 


(nt)’ 


P n ( 0 ; 0 , 7 ) “‘ x ' ' 2 “ v ”' ' 4 

in accordance to ( 13 . 731 ) . From these facts c 3 is identify to be 

hi (ro + 7_A-y)) _ i)i„ 7 + 7 _ii n(2 7). 

Consequently, the equation (j 3 . 74 [) becomes, 
c 2 (q + 1 ) - c 2 (7) = c 3 (7) = In ^ 1 + 7 ^ A ——^ - ^ 7-0 In7 - 7 + i In (2 tt) 


We may identify 

r-„l+A = In f 

r(7)G( 7 ) 

up to another constant b( 7), which satisfies, 


, , , (G( I + 7-A-/3)'. ,, . 

c 2 (7) = In - — 7777 . -- + 6(7), 7 > 0 


6(7 +1) - 6(7) = f ™ ( 1 -- + 


2 t e* — 1 It 


=-7 1 


-dt. 
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Binet’s formula, see ([39], P249), says that, 

1,1 rW = ( z - 1 ) lll(z) - z + \ ln (2l,) +1 G - T + ?=l) V*’ 

Remark 8: As 7 —> 00 , 


Re z > 0. 


6(7 + 1) - K7) = ^ - 3^5 + 0(l/7 5 ). 


We also note that, 

6(7) = «-^ln 7 +^-^ + +°(l/7 4 ). ( 3 . 80 ) 

It may be of interest to determine the constant k . 


4 Limiting behavior of the Kernel with the Pollaczek- 
Jacobi type weight. 


The ladder operator formalism for orthogonal polynomials and the associated compatibility 
condition can be found, for example, in [IT] . This is a handy tool regarding orthogonal poly¬ 
nomials on the unit circle [3j, construction of Jacobi polynomials mi and has been adapted 
to obtain the Painleve equations, arising the deformation of classical weight. We refer the 
Reader to 0 HOI U m\ [ 221 12 H]. In particular, this r formalism has been applied to spe¬ 
cific orthogonal polynomial ensemble in [HI EE], and its comparison with the isomonodromy 
approachp?]. For the extension of ladder operators to discrete orthogonal polynomials, and 
allied discrete Painleve equations; to q -orthogonal polynomials and the allied q -Painleve 
equations, we refer the Reader to [ 6 , 8 j [TJ [16] i3U [32] 12]. 

In this section, we construct the Lax pair for our problem. 

First, we recall the ladder operator relations of (2.7) and (2.8) in |10j . satisfied by the 
polynomials P n (x]t,a, (3) orthogonal with respect to the Plooaczek-Jacobi type weight, 


+ B n (x ; t, a , /?) ) P n (x ; t , a, (3) = /3 n A n (x ; t, a, /3)P„_i(x; t, a, /3), 


(4.81) 


(T - B n (x\t,a,fi) - v'(i)) P„-i(x-,t,a,p) = -A n _i(x\t,a,p)P n (x\t,a,P), (4.82) 

The “coefficients”, A n (x]t,a, (3) and B n {x\ t, a, (3) are given by (2.14), (2.15) in [TO] , 
respectively; 


. , R* R n 

A n (x;t,a,/3 ) = — + 


x* 


X x — 1 X z X x — 1 


and R* n , R n , r* , r n are defined by (2.16) - (2.19) in [TO], all depending on t, a, /3. 
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Theorem 12. The monic polynomials P n (x\ t, a, (3) orthogonal with respect to the Pollaczek- 
Jacohi type weight, satisfy the ladder operator relations in t, 


xt— - r* ) P n (x ; t, a, (3) = /3 n R* n P n _ i(x; t, a , (3), 


xt^ + r* - t + xR* n _^j P n ^i(x-,t,a, f3) = R^P^x^a, f3). 
By the Christoffel-Darboux formula [43], the reproducing kernel reads, 


(4.83) 

(4.84) 


K n (x,y) := 


K 0 n (x)(/) n _i(r/) - 0 n (y)0 n _i(a;) 
h n —i x-y 


(4.85) 


where h n is the square of the L 2 norm is given by (II.3p and (3 n = h n /h n _ Here, 4> n (x) 
and 0 n _i(x) are defined in terms of P n (x]t,a, (3) by 

, P n (x\t,a,P ) i, . P n (x\t,a,P) <*, .4 t 

(x) -;=- wz(x-,t,a, p) — - -= -# 2(1 — x) 2 e ^, 


.-i W : = 


\/hn 

P n -i(x', t, a, (3) 1 


y/hn 


\j h n — i 


Let 4/ , 


vh = 


X 




, r,\ Pn—l( x R]Qtil3') 2 ./.. \ S. -X- 

W2(x]t,a, p) — - . - £2(1— £) 2 e a*. 

yhn -1 


P n (x:t,a,/3) —/i \£ —£- 

X 2 (1 — XJ 2 e 

P n —i(x:t,a,8) —/-< —— 

h-— —£2 1 — x be 2 b. 


We find, 


and 


.Ci a c 3 . T 

4b = I -- + — + -4 ) 4h 

£ — 1 £ £ z 


Ci = 


h t 0 \/ 3n /~i 

- , h'2 — 


Rn—1 \J ftn 8 n 


c 3 = 


2 J 


n r n 4 j Rn\/(3n 

Rn—l^ZPn ^ 4“ y 


-K-iVfa -| + < 


(4.86) 


(4.87) 


(4.88) 


Similarly, the ladder operator relationships in t, see (j4.83j) and (j4.84j) can be re-written as 


and 


4b = (A) + tf, 


Rn 

1 

7 

O ttf 


r < i 

RnV n 

2 t 

0 

, D 1 = 

t 2 

K-iVdh 

t 

1 < 

2 1 


t 

2 t . 


(4.89) 


(4.90) 
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The equations (j4.86j) and (I4.89j> give the Lax pair mentioned earlier. The compatibility 
condition of the Lax pair reads = ^> tx , which, when un-pact, re-produces the differential 
differences obtained in ma- 

Writing out (14.861) in component form, followed by eliminating 0 n _i(x) in favor of (f> n {x) 1 
we find, 


2 ( x 3 — x 2 ) 0"(x) + 2 ( 3x 2 — 2x ) 4>' n {x) — (2(2 n + a + f3)x + 2 (r n — r*) + t — 2n — a) 0^ 

Rt - Rn 


2 (R n Rn) x ( X jO -/ / \ , 

~ ^ ~ K ^ + 


((2 n + a + j3)x 2 + (2(r n - r*) 


x(i?* - - R 

+t-2n-a)x + 2r* n -t)<f> n - 2x 2 (l_i\ [( x ( 2 ( r ^ “ r n) + t - 2n - a) + 2r* - t 
+x 2 {2n + a + /3)) 2 - 4 (z(.R* - -R„) - i?*) (x(i2*_i - R n -i) ~ R* n - 1 ) Pn 4>n = 0, (4.91) 


where r n , r* , R* n , i? n , i?*_ x , i? n _i may be expressed by H n , H' n (t ) and H n ”(t), these 
are (2.56), (2.32), (2.33), (3.7), (5.1), (5.2), (5.8) in [TO] . For convenience, we list these 
relations in Appendix C. 

We give in the next theorem, the scaled kernel. 

Theorem 13. Let 

C:=4n 2 x, C*:=4n 2 r/, s\=2n 2 t, 0(C) := lim 0 n (4.92) 

MOO 04?!^/ 


where n —* oo, t —> 0 + , C an d C* compact subsets of (0, oo), and s finite, then, 


lim —K n 

imoo 2n \ 4n 2 4n 2 / 

(C - 4aW')0(C)C*V(C*) - (C* - 4.s77')0(C*)Cy (C) - 2 S 2 77"(C - C*)0(C)0(C*) M v 

(C-C*)(C-4aW)(C*-W) ’ J 


Furthermore, 0(C) satisfies the ODE, 


CV'(C) + (c 


4att , QQC \ 
C-4 S ?7'( S )j 


¥0 + 


C 2s 2 77"(s) 

4 C - 4sW(s) 


4 C 2 


0(0 = 0 , 

(4.94) 


with the boundary conditions of 0(0, ex, 0) = 0, a > 0 . 

Proof. 0 n _i(x) may be expressed in terms of 0^(x) and 0 n (a;) via one of the ladder operator 
relations, 

, . _ ((2n + a + 0)x 2 + (2(n + r n - r*) + t - a)x + 2r* - t) 0 n - 2(x - l)x 2 0 / ri (a;) 

“ 2.yJ~fif{xR n + (1 — %)R*i) ' 
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Substituting 0 n _ x , given above, (14.92 p and the large n asymptotic of r* , r n , R n and R* n 
presented in the Appendix C, into the kernel given by the equation (I4.85f) . and let n —* oo, 
then the limiting kernel (j4.93j) is obtained. 

We first substitute ( = 4 n 2 x and s = 2 n 2 t into (14.911) . satisfied by 0 n (:r). Recall r n , 
r* , R* n , R, n , R n —i in [TO], listed in the Appendix C. Now let n —> oo, then 0(C) 

satisfies the second order ODE, namely, (14.941) . □ 

We list here some properties of the limiting kernel. 

Let s = 0. The boundary condition, Q,a,f3) = 0, implies the kernel reduces to 


lim —K„ 

iho o 2 n 


4n 2 ’ An 2 ) 


meeto - ckoc^o 

c-c* 


In which case, the equation (14.941) becomes the Bessel differential equation, with the bound¬ 
ary condition of 0(0) = 0. Then the limiting kernel becomes the celebrated Bessel Kernel 
of Tracy and Widom m ■ 

For s > 0, we make the transformation, 


C^4sH'( s )C, C* -> 4sM'(s)C*. 0(0 -> C^'VK). 
then the kernel (14.931) becomes as, 


1 f C C 

lim —IfC, ' 


A(C)B(C)-A(e)B(0 


>oo 2n n \An 2 'An 2 1 
where the functions A(C) and B(() are given by 


C-C* 


fp(s) £ P (s) (2 

A( C) := — 4 === 0(C), 5(C):=-4==7^-r0'(C), 


2,/7H 7 {s) 


2^/7W{s) C-l 


and 


' ,(S) := 2Ws) 


The ODE (14.941) transforms to, 

2(1 + bo) 1 


</>”(() + 


C C-l 


)^c) + (-b-j| + L + h) 0( C )-o, 


where 


a 


{sW'f - - a 2 W 2 

AW 


(4.95) 


b 4 = —(nr, b 3 = -, b 2 = 

16 V ; ’ AW 


sW 


-, 6, = sft', 5 0 = ^7- 


The coefficients 6 X , b 2 , 63 , 64 satisfy the relation, 


61 + &2 - 63 — h = 0, 
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which is the a-form of the Painleve equation found in the Theorem 3. 


(. sH"f + 4 {Wf {sW - H) - (aW + 0 =0. 

The kernel (14.951) has continuity property, a feature of the integrable kernels, since it can 
be written as, 


lim — 

n—>oo 2 n 


c e 


4 n 2 ’ 4n 2 


E}=o/i(0 9j(C) 
C-C* 


where 


sU" 

0 . r w 

MO = (-i) 2 " J — 


2y / kH 7 VC-l^C 


C 2 dV (w 

^ ' 0(0, <7,(0= - 


C 2 d 


2 OH 7 VC - 


1 -,' 


0 ( 0 , 


and j G {0,1}, with the property ^]=o MOM 0 = 0 . 

Using the continuity property, the limiting kernel (j4.93[) becomes, 


lim — 

n—>oo 2n 


c c 


4n 2 ’ 4n 2 


C 2 0 / (O - 00(0000 - 2(0(00(0 00(00(0 - 2s 2 n"(s)^(0 

C-4 sH'(s) (C-4sn'(s)) 2 

c 2 ( 0 (c)) 2 + (5-T-p-?- w ) ^(0 


C 2 (0'(O) 2 + 


c - 4sW 

J - £ ) 0 2 (O 


s > 0, 


c 


s = 0, 


= A 


J'o(VC) - Jo .-1 (VC) J a+l (Vd 


where J a (z) is the Bessel function of order a , A is a parameter. The first equality above, is 
found by applying L’Hospital rule to (14.93 1) , the second equality is achieved by eliminating 
0"(C) with the aid of (14.941) and "H(0,a,/3) = 0. The Bessel equation satisfied by 0(C) with 
the boundary condition 0(0) = 0 has regular solution 0(C) = A J a {VC ) with a parameter 
of A. 


5 Limiting Kernel with the perturbed Laguerre weight. 

In this section, we make use of the method in the previous section, and apply to the kernel 
generated by the singularly perturbed Laguerre weight. It is interesting that the limiting 
kernel via from the perturbed Laguerre weight is the same with the limiting kernel that arises 
from the Pollaczek Jacobi type weight, although their scaling schemes are quite different. 
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First, the monic orthogonal polynomials P n {x\t,a) with respect to the singular per¬ 
turbed Laguerre weight, 

w(x;t,a) = x a e~ x e~^, 0 < x < oo, a > 0, t> 0, 

satisfy the ladder operator relations have been derived in |18j , which we restate here, 


^- + B n {x) ) P n (x;t,a) = f3 n A n (x)P n _i(x] t, a), 


(jL ~ Bn ^ Pn-i(x;t,a) = —A n - 1 (x)P n (x; t, a), 

where v(x) = — In w(x; t, a ). The coefficients A n (x), B n (x) are given by (2.7) and (2.8) 


m c _ a 

A ( \ _ 1 . Rn r> / \ _ n | r n 

A n [x) — —I-y, B n {x) —- 1 — 5 -. 

We use Rn(t), r n (t) and t instead of a n (s ), b n (s) and s respectively in [T8] ; these are 
defined by, 


t 


Rn(t):=— -P n (y;t,a)w(y;t,a)dy , r n (t) : = 

"n Jo V 


t 


hn— 1 Jo V 


Pn(y; t, a)P n -i(y; t, a)w(y)dy, 


depend on t and a, with the initial conditions R n ( 0) = 0, r n (0) = 0. 

The monic orthogonal polynomials P n (x]t,a ) satisfy ladder operator relations in t, see 
(5.56) and (5.57) in [18], 

d Tn \ / 4 . \ finRn D / , \ 

Jt~Vt) Pn ^ X] t} a > = - B~ Pn ~^ X] a 

d 1 r n R n -i \ „ , . R n ~ 1 „ , , 

—-1—- H-— P n - i(aj; t, «) = —— P n (x ; t, a). 

dt x xt t J xt 

The reproducing Kernel with respect to w(x;t,a ) is given by, 

1 

K \ 2 tp n {x)ip n -i{y) - <p n {y)(p n - i(x) 


K n(x,y) := 


hn —1 


x-y 


(5.96) 


where, 


/ \ Pn(^£} t , Cl) 1/ \ PnijE 1 11 OtJ) a _x t_ 

(p n [x) := - -f= — w 2 (x) = - n —x 2 e 2 2 ^, 


yfhj 


y/hj 


/ \ Pn—lijE , t, Cl) i» , P n —\{x^t ) odj a _ x _t_ 

<^ n _i(a;) :=-— w*{x) = - -== — x*e 2 2 -, 

yh n -1 \/h n -1 

where h n is the square of the L 2 norm and f3 n = h n /h n -i ■ 
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Let <f>(x) 


<f> := 


<7>n(x) 


<£ n l(x) 



We find, 

where dq denotes , 

a, r;rr 7 

1 


Pn(x:t,a) — — £—L ' 

nV J-i 2 e 2 2* 

V Rn 

P n -i(x:t,a) SL _£__L 
i-J-L 2 e 2 2 x 

‘sj ^Xn — 1 


= ( £o + — + ^ 
x x z 


(5.97) 


£n = 


-1 0 

0 1 


E 1 = 


n+ § ^ 


The quantity $ also satisfies, 


where 


\ffin 

~ n ~ fj 

$t = 

R + h 

V x 


el = 


2 Rn \J (3 n 

Rn— lVS 2 ^ 


(5.98) 


^0 = 

0 

2t ;; 

m Rn — 1 

, E\ = 

la 1 

t 2 

Pn \J fin 

t 

1 7'n 


L U 2i J 


t 

2 t J 


The equations (15.971) and (I5.98P is the Lax pair of this problem. Un-packing the compatibility 
condition is satisfied by the Lax pair, we find the following set of scalar equations, 


d 2 (r 2 — tr n ) , . 

t~n r n = -5-- + (2n + 1 + ct)r n - nt , 

Chi' -LLyi 

t—R n = R 2 + (2n + 1 + a)R n — t + 2r n , 
t-^_R n -1 = t - 2r n - (2n - 1 + a + R n _i).R n _i. 


(5.99) 

(5.100) 


Note that the above two Riccati equations are the same with (3.10) and (3.11) in [IS] , the 
third equation in the above also can be derived by (2.16) and (3.10) in [18j. 

We write (15.971) as a set of scalar equations, 


1 2 n + a t — 2r n \ . . (1 R n , ^ 

Vn( X ) = ( ~2 + - 3 ^“ + + ( - + ) PZVn-l(x), 


<Pn -lW = 


1 R„_i\ | . . / 1 2n + a f-2r n 

H-— ) Pn^n{x) ~ -- H---h 


X 


X z 


2x 


2x 2 


If we eliminate tp n _i(x) The resulting ODE reads, 

x — R 


2 n/ \ , / i xR n 

X <p n (x) + [x + 


X + 

2 t 2 




fd_ TD _ 77> 

n 1 dt 1X,n 1X,n 

2(x + R n ) 


Vn-l(x). 


<Pn(x) 


(5.101) 

(5.102) 


x" a" t“ ( a\ t at „ . . . 

T + T + 4^“ I (" + 1+ 2)“2 + i; + - ff d V’nW = 


(5.103) 
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with the aid of 0 5.100 1) , and (2.12), (2.13), (3.21), (3.29) in 
In order to investigate the limiting behavior of the kernel, we recall the equations of (2.12), 
(2.13), (3.21), (3.29) in 


The equations, (3.28) in |il8j becomes, 

oi n — 2n T 1 T ol T 
and with (2.9) in [T8], we find 

Rn = 

Let 


/3 n (R n + R n - i ) = -(2 n + a)r n + nt, 
finRnRn—l ^n ^n ■ 

I3 n = n 2 + na + tH' n - H n . 
2 - tH ' n ) 


T — / 

'n L11 ni 


tH + n — (2 n + a)H' n ’ 


2 (t(^; ) 2 - ^;) 


tH'[ + n — (2n + a)H' n ’ 
"H(s,a) := lim 


(5.104) 


a 


2n + 1 + a 

In the next theorem we describe the scaled kernel. 

Theorem 14. Let 


£ := 4na;, £* := 4 ny, s := (2n + 1 + a)t, and </?(£) := lim <p n . , , , 

n—> oo \ 472 / 


(5.105) 


where t —)• 0 + , n —>■ oo and £, £* are in compact subsets of (0, oo), s finite, a > 0, then 

lim —K n f— 

n—too 4n \4n 4n 

(j - 4s?^)M0ey(r) - (r - 4sK'(s)Moey (Q - 2s 2 ?r(s)(e - rMeMO 

(e-e)(e- w( S ))(£*-4sw(5)) 

and <^(£) satisfies the second order ODE, 

4s77'(s)£ \ , /£ , 2s 2 ?7"(s) 


(5.106) 


ev(o + e 


£ - 4sH'{s) 




a 2 s 2 


as 


£-4sft'(s) 4 £ 2 


— - n(s) \ m= 0, 

(5.107) 


with the boundary conditions 99 ( 0 ) = 0, 
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Proof. By (15.1011) . ip n _i(x) can be represented by <p n (x) and its derivative as, 


<Pn-i(x) 





Inserting the above equation and (15. 105ft into the kernel given by (15.96ft , let n —> oo, then 
we arrive at (15.1061) and (15.1071) . 

In deriving these results, a number asymptotic relations, are listed here. 


r n 


tH' n (t) ~ sPL'(s), 


Rn 


r\j 


sPL'(s) 

n 


t dt Rr ‘ 


■sPL'(s) 

n 


s 2 PL"(s) 

n 


where t —^ 0 , n —> oo and s = ( 2 n + 1 + a)t is hxed. 


□ 


Hence, we see that the scaling limit of the logarithmic derivative of the Hankel determi¬ 
nant H n (t, a) generated by the perturbed Laguerre weight is the same with the scaling limit 
of logarithmic derivative of the Hankel determinant generated by the Pollaczek-Jacobi type 
kernel. Moreover, the limiting kernel (1 5.106 ]) is completely characterized by the second 
order linear ODE of (1 5.1071) . which is the same with the limiting Pollaczek-Jacobi type 
kernel in the Theorem 13, both of which are “scaled at the origin” but their scaling scheme 
are different from each other. 

We present here another version of the ODE (1 5.107 1) . 

For s > 0, we make the transformation, 


£ 4 sH% 


where PC denotes dPL(s)/ds , into (1 5.107 1) : we see that </?(£) satisfies, 

p"K)+ (I - ^ T ) m +(-| -1+p+ f +-^u) rn= o. 


where, 


a 


a 4 = 


16^/2’ “3 •- 4 2 2 PC 

The coefficients of a 4 , a 3 , a 2 , aq and a 0 satishes the relation, 

<3q + 2a 0 + cq + CI2 — — 04 = 0, 

which is the a -form of the Painlcve equation, 


a 2 = —PL 


a 

T 


sW‘ 


a 4 = sPL' 


sPi h 


(sPi"f + 4 (PL') 2 (sPC - PL) - (olW + 


= 0. 


sPi" 


a 0 — 


2 PC' 


We consider an interesting special case, where a = 0 The ODE becomes 


</(£) + ( j 


£ + 2 si 


¥>'(£) + -71 + 


s“ 27ss + 5 


£ 4 


36£ 2 


+ 


3 


+ 


12ss£ 12(^ S 3+2 s3) 


¥>(£) = 0 , 
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with the boundary condition of <p(0) = 0. The solution of the a-from Painleve equation, in 
this situation, reads, 


„„ , 3 2 1 

WW--j*» + 36- 

The region £ << 1, and s » 1. The approximating equation is 

V’”K) + |v>'K)-^tt) = o, 

with boundary condition -0(0) = 0. The solution, up to a constant multiplier, is given by 

_s 

e ?. 


The region £ >> l,£s s >> 1. The approximating equation is 


no + -^no + 


3ss — 1 

— —I - 


12ss£ 12 ss(£ + 2s3) 


m = o, 


we find that, <£>(£) is asymptotic to 

d (£ + 2s§)HeunC ^0,1, 0, + 1^ , 

where A is an arbitrary constant. 
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Appendix A 

We list here a selection of integral formulas that are relevant for the computations in the 
main text. These can also be found in [20J and |T3]. For 0 < a < b, we have, 



y/W~ x )( x - a ) 


dx = 7r. 


x 


yf{b — x) (x — a) 


dx = 


(a + b)7T 



Xy/{b 


V(J> 


==dx 

— x)(x — a) 


= dx 

x)(x — a) 


7r 

Vab 

(a + b) 7T 

2(a6)§ 


(Al) 

(A2) 

(A3) 

(A4) 
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y/a + \fb ^ 

(A5) 

y/(b — x)(x — 

a) \ 

, 2 J 

Inx 

dx- 2n 

2 y/ab 

(-46) 

X\J (b — x)(x 

— a) v ab 

V^+Vb' 


Appendix B 

We list here a number of identities involving R n , R n *, r n , r n *,(3 n , H n which can be 
found in [TO] , 

R* n = R n — (271 + 1 + ol + /?), 


Kf ~ K = p n R* n R* n _ A 
f 'n T / f f // Rn Rn — I ■ 


-« - r n ) 2 - (/3 + t)r n + (f - a - 2n)r* + nt 
1 — (2n + a + /3) 2 

* = nt + tH' n 
n 2n + a + /3 ’ 


= 


n(n + a) + tH' n — H n 


Rn = 


2n + a + f3 

(2 n + 1 + a + /3) [2r 2 + (f + 2/3 - 2r* n )r n + (2n + a)r* 


nt - tr' n ] 


2[(r* - r n ) 2 + (2n + a - t)r* + (/3 + t)r n - nt] 
The expressions of these for large n , are given below, 


<(*) 


sH! s — (a + /3)sK' 
2 n n 2 



n a — f3 f5 2 — a 2 — A(T-i — sW) (a +/3)((a 2 —/3 2 ) + 4 (Ji — sH)) 1 

’■"=2 + — + -ST-^ + i - 16tf -^ + 0( ^' 

n n z n 6 


2sW (a + f3 — l)sW - 2s 2 U" ^ 1 , 

K =-+ -—-—^-+ O(-), 


n 


n* 


n° 
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where PL is given by 


PL(s,a,/3) := lim H n (-^-,a,/3) , 
n-¥ oo \2n z / 

PL' denotes dH(s)/ds and s = 2n 2 t. 
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